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Abstract. In the present paper, we prove the a priori estimates of Sobolev norms 
for a free boundary problem of the incompressible inviscid MHD equations in all 
physical spatial dimensions n — 2 and 3 by adopting a geometrical point of view 
used in [Q , and estimating quantities such as the second fundamental form and the 
| velocity of the free surface. We identify the well-posedness condition that the outer 

■ normal derivative of the total pressure including the fluid and magnetic pressures 



is negative on the free boundary, which is similar to the physical condition (Taylor 



^ sign condition) for the incompressible Euler equations of fluids. 
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3 '. 1- Introduction 

m ■ 

In the present paper, we consider the following incompressible inviscid magneto- 
hydrodynamics (MHD) equations 

v t + vdv + dp = — (B-dB- -d\B\ 2 ), in 0, (1.1a) 

*jj ' 4:71 2 

B t + vdB = B-dv, in & (1.1b) 
divw = 0, div5 = 0, (1.1c) 

describing the motion of conducting fluids in an electromagnetic field, where the 
velocity field of the fluids v = (v i, • • • , v n ) , the magnetic field B = (B\, • • • , B n ) 
and the domain @ C [0, T] x R n are the unknowns to be determined. Here n e {2, 3} 
is the spatial dimension, p is the fluid pressure, d = (di, • • • ,d n ) and div are the 
usual gradient operator and spatial divergence. Given a simply connected bounded 
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domain Sl^ C W 1 and the initial data Vq and B$ satisfying the constraints div v o = 
and div Bq = 0, we want to find a set 3i C [0, T] x M n and the vector fields v and .B 
solving ( |1 . 1| ) and satisfying the initial conditions: 

{x : (0, x) G £>} = %, {v, B) = {v , B ) on {0} x %. (1.2) 

Let S) t = {x E W 1 : (t,x) G we also require the following boundary conditions 
on the free boundary d@ t : 

Vj\f = k on d@ t , (1.3a) 

p = ond@ t , (1.3b) 

\B\=q and B-N = on d@ u (l-3c) 

for each £ G [0,T], where A/" is the exterior unit normal to d£t t , = ^™ =1 A/\>i, 
and k is the normal velocity of d&ti ? is a non-negative constant. Condition (|1.3c|) 
should be understood as the constraints on the initial data. Indeed, we will verify 
that the condition B ■ N = on d% holds for all t G [0, T] if it holds initially. We 
remark here some physical meaning of the boundary conditions. Condition ( |1.3a|) 



means that the boundary of % moves with the fluids, ( |Ob|) means that outside 



the fluid region @ t is the vacuum, the condition B ■ N — on d^ t comes from 
the assumption that the boundary dS> is a perfect conductor. Indeed, if we use 
E to denote the electric field induced by the magnetic field B, then the boundary 
condition B-M = on dSt t gives rise to E x M = on d@ t . The boundary condition 
\B\ = const on d£t t (the magnetic strength is constant on the boundary) is needed 
to guarantee that the total energy of the system is conserved, i.e., 

d f fl, ,o 1 



dt 



jf \^\v\ 2 + ^\B\ 2 j (t,x)dx = 0. 



Condition (|1.3c[) includes the widely used (e.g., |12[) zero magnetic field boundary 



condition as the special case, but it is much more general and physically reasonable. 

In the classical plasma- vacuum interface problem (cf. [^,^3[), suppose that the 
interface between the plasma region Q p (t) and the vacuum region Q v (t) is T(t) which 
moves with the plasma, then it requires that ( |1 . 1| ) holds in the plasma region Q p (t), 
while in the vacuum region Q v (t), the vacuum magnetic field SS satisfies 

VxJ = 0, V-^ = 0. (1.4) 

On the interface T(t), it holds that 

p = 0, \B\ = \&\, B-N ' = SS-N ' = 0, (1.5) 

where is the unit normal to T(t). Therefore, the boundary conditions in ( p..3| ) 
also model the plasma- vacuum problem for the case when constant. 

We will prove a priori bounds for the free boundary problem (|1 . 1|) , ( |1.2|) and (1.3) 
in Sobolev spaces under the following condition 

V ^ (p + ^l 5 ! 2 ) ^- £<0on dS) t, (1.6) 

where V^f = M l di. We assume that this condition holds initially, and will verify 
that it holds true for some time. For the free boundary problem of motion of incom- 
pressible fluids in vacuum, without magnetic fields, the natural physical condition 

(cf. i,il,|,|iM3M^llll ; 0) reads that 

^MV < -e < on d& t , (1.7) 
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which excludes the possibility of the Rayleigh- Taylor type instability (see ||). In 



this paper, we find that the natural physical condition is (|1.6|) when the equations of 
magnetic field couple with the fluids equation. In fact, the quantity p+ -^\B\ 2 , the 
total pressure of the system, will play an important role in our analysis. Roughly 
speaking, the velocity tells the boundary where to move, and the boundary is the 
level set of the total pressure that determines the acceleration. 

The free surface problem of the incompressible Euler equations of fluids has at- 
tracted much attention in the recent decades. Important progress has been made 
for flows with or without vorticity, with or without surface tension. We refer readers 
to [§|,|,|,|T§^gl|,g4|,0,g3. 



On the other hand, there have been only few results on the interface problems for 
the MHD equations. This is due to the difficulties caused by the strong coupling 
between the velocity fields and magnetic fields. In this direction, the well-posedness 
of a linearized compressible plasma-vacuum interface problem was investigated in 



23[ , and a stationary problem was studied in || . The current- vortex sheets problem 
was studied in || and [p2| . For the incompressible viscous MHD equations, a free 
boundary problem in a simply connected domain of M 3 was studied by a linearization 



technique and the construction of a sequence of successive approximations in [18 
with an irrotational condition for magnetic fields in a part of the domain. 



In this paper, we prove the a priori estimates for the free boundary problem ( 1.1 
( |1.2|) and ( |1.3| ) in all physical spatial dimensions n = 2, 3 by adopting a geometrical 
point of view used in ||, and estimating quantities such as the second fundamental 
form and the velocity of the free surface. Throughout the paper, we use the Einstein 
summation convention, that is, when an index variable appears twice in a single 
term it implies summation of that term over all the values of the index. Denote the 
material derivative D t = d t + v ■ d and the total pressure P = p + g^|I?| 2 , we can 
write the free boundary problem as 

D t Vj + djP = —B k d k B j in 9, (1.8a) 
An 

D t Bj = B k d kVj in 9, (1.8b) 

djV j = in 9\ djB j = on {t = 0} x (1.8c) 

v n = k on [0,T] x d%, (1.8d) 

\B\ = q on d@, BjM j = on {t = 0} x d%, (1.8e) 

p = on ^, (1.8f) 

VmP < on {t = 0} x 8%. (1.8g) 

We will derive the energy estimates from which the Sobolev norms of H s (£$ t ) (s ^ 
n + 1) of solutions will be derived. For this purpose, we define the energy norms as 
follows: The zeroth-order energy, E (t), is defined as the total energy of the system, 
i.e., 

E {t)= f 6 ij (v i v j + -l-B i B j )dx, (1.9) 

which is conserved, i.e., 

E (t) = E (0), for ^ t ^ T. (1.10) 

The higher order energy norm has a boundary part and an interior part. The bound- 
ary part controls the norms of the second fundamental form of the free surface, the 



4 



CHENGCHUN HAO AND T. LUO 



interior part controls the norms of the velocity, magnetic fields and hence the pres- 
sure. We will prove that the time derivatives of the energy norms are controlled by 
themselves. A crucial point in the construction of the higher order energy norms 
is that the time derivatives of the interior parts will, after integrating by parts, 
contribute some boundary terms that cancel the leading-order terms in the corre- 
sponding time derivatives of the boundary integrals. To this end, we need to project 
the equations for the total pressure P = p + g^r|-B| 2 to the tangent space of the 
boundary. The orthogonal projection IT to the tangent space of the boundary of a 
(0, r) tensor a is defined to be the projection of each component along the normal: 

(Ua) iv .. ir = nj • ••Ii;;n. r . ,,. where U{ = 5{ - (1.11) 

with J\f j = 5 ij Mi = Mj. 

Let di = U^dj be a tangential derivative. If q = const on d%, it follows that 
diq = there and 

ilid 2 q) l3 = 6 tJ V X q, (1.12) 

where 0y = dj\fj is the second fundamental form of dS> t - The higher order energies 
are defined as: For r > 1 

E T (t) = J 5 ij (Q(d r Vi ,d r Vj ) + ^-Q(d r B i ,d r B j )J dx 

+ J (V^curltf + ^ld^curlEl 2 ^ dx 

+ I(r) [ Q(d r P,d r P)$dS, (1.13) 

Jd&t 

where J(r) = if r = 1 and I(r) = 1 for r > 1, so we do not need the boundary 
integral for r = 1, 

■& = (-VatP)- 1 . 

Here Q is a positive definite quadratic form which, when restricted to the boundary, 
is the inner product of the tangential components Q(at,(3) = (Ha, 11/3) and in the 
interior Q(a,a) increases to the norm \a\ 2 . To be more specific, let 

Q(a, 0) = q^ ■ ■ ■ q trjr a n ... lr (3 h ... jr (1.14) 

where 

tfi = #3 - n (dfM i U j , d(x) = dist (x, d%), M i = -5 ij djd. (1.15) 

Here n is a smooth cutoff function satisfying ^ n(d) ^ 1, n(d) = 1 when d < do/4 
and n(d) = when d > do/2, do is a fixed number that is smaller than the injectivity 
radius of the normal exponential map i , defined to be the largest number i such 
that the map 

d% x (-t 0j t ) -} {i G 1" : dist (x, d%) < l } (1.16) 

given by 

(x, i) — y x — x + lM{x) 

is an injection. 

The main theorems in this paper are as follows: 

Theorem 1.1. For any smooth solution of the free boundary problem (|1.8|) for 
^ t ^ T satisfying 

\dP\ ^ M, \dv\ <M, in @ u (1.17) 
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|0| + \dv\ + - ^K, on d%, (1.18) 

we have for t G [0, T] 

#i(*) < 2e CMt £i(0) + CK 2 (Vol ^ + E (0)) [e cm - l) , (1.19) 
for some positive constants C and M . 

Theorem 1.2. Let r G {2, ••• ,n + 1}, t/ien t/iere exists a T > suc/i i/ia£ £/ie 
following holds: For any smooth solution of the free boundary problem ( |1.8| ) /or 
^ t ^ T satisfying 



\B\ ^Mi 


for r = 2, 


in 




(1.20) 


|<9P| < M, <M, 


\dB\ ^ M, 


in 




(1.21) 


|0| + 1/io 




on 


d@ u 


(1.22) 


-VatP ^ £ >0, 




on 


d%, 


(1.23) 


\d 2 P\ + \V*D t P\ ^L, 




on 




(1.24) 



we have, for t G [0, T], 

E r (t) ^ e Clt E r (0) + C 2 (e Clt - l) , (1.25) 

where the positive constants C\ and C 2 depend on K, K 1 , M, M 1; L, 1/e, Yo\@> t , 
E o (0), E^O), andE r ^(0). 

Most of the a priori bounds (|1.20|) - (|1.24|) can be obtained from the energy norms 
by the elliptic estimates which are used to control all components of d r v , d r B and 
d r p from the tangential components Ud r P in the energy norms, and a bound for 
the second fundamental form of the free boundary 

\\B r - 2 e\\ L 2 im) < C f#,L,M, J,^, Vol W) E r 

for r > 2, which controls the regularity of the free boundary. 

Since E (t) = E (0) and Vol % = Vol^o, recursively we can prove the following 



main theorem from Theorems |1.1| - |L2 . 
Theorem 1.3. Let 

/C(0) = max (||0(O, .)|| L -. (Mb) , V*o(0)) , (1-26) 

5(0) = || 1/(V* P(0, ■)) || Loo(a3b) = 1/6(0) > 0. (1.27) 

There exists a continuous function T > such that if 

T < T(/C(0), 5(0), E {0), ■■■ , P„ +1 (0), Vol %), (1.28) 

then any smooth solution of the free boundary problem for MHD equations (|1.8j ) for 
^ t ^ T satisfies 

n+l n+1 

E s (t) < 2 E s (0), O^t^T. (1.29) 

In order to prove the above theorems, we need to use the elliptic estimates of the 
pressure p. However, the time derivative of Ap involves a third-order term of the 
velocity which needs to be controlled by higher order energies. In order to overcome 
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this difficulty, we work on the equations for the total pressure P = p + ^|5| 2 , 
instead of those for the fluid pressure p. 

Before we close this introduction, we mention here some studies on viscous or 
inviscid MHD equations, including the Cauchy problem or initial boundary value 
problems for the fixed boundaries P, |TTHT4l, P^8|-pO|, [26|j and the references therein. 

The rest of this paper is organized as follows: In section 2, we use the Lagrangian 
coordinates to transform the free boundary problem to a fixed initial boundary prob- 
lem. The Lagrangian transformation induces a Riemannian metric on 3t , for which 
we recall the time evolution properties derived in H, and prove some new identities 
which will be used later. We also write the equations in Lagrangian coordinates, 
by using the covariant spatial derivatives with respect to the Riemannian metric in- 
duced by the Lagrangian transformation, instead of using the ordinary derivatives. 
In section 3, we prove the conservation of the zeroth order energy Eo(t), from which 
one can see that the boundary conditions on the magnetic fields B is necessary for 
this energy conservation. We also prove in section 3 that the condition B • J\f = 
on the boundary propagates along the boundary. Section 4 is devoted to the first 
order energy estimates. In section 5, we prove the higher order energy estimates, by 
using the identities derived in section 2, the time evolution property of the metric 
on the boundary induced by the above mentioned Riemannian metric induced by 
the Lagrangian transformation, the projection properties and the elliptic estimates. 
In the derivation of the higher order energy estimates in section 5, some a priori 
assumptions are made, which will be justified in section 6. We also give an appendix 
on some estimates used in the previous sections, which are basically proved in [Q. 

2. Reformulation in Lagrangian Coordinates 

Assume that we are given a velocity vector field v(t,x) defined in a set ^ C 
[0, T] x R n such that the boundary of % = {x : (t, x) G 2>} moves with the 
velocity, i.e., (l,v) G T{d$>). We will now introduce Lagrangian or co-moving 
coordinates, that is, coordinates that are constant along the integral curves of the 
velocity vector field so that the boundary becomes fixed in these coordinates (cf. 0]). 
Let x = x(t, y) = ft(y) be the trajectory of the fluid given by 

dx 

— = v(t,x(t,y)), (t, y) G [0, T] x Q, 

x(0,y) = fo(y), y e n. 

where, when t = 0, we can start with either the Euclidean coordinates in Q = £? 
or some other coordinates /o : Q — > where /o is a diffeomorphism in which the 
domain Q becomes simple. For each t, we will then have a change of coordinates 
ft : ft — > @u taking y — > x(t,y). The Euclidean metric 5ij in @ t then induces a 
metric 



Qx % Ox^ 

9a^y) = 5 lJ —— b (2.2) 



and its inverse 



in Q for each fixed t. 

We will use covariant differentiation in Q with respect to the metric g a b(t, y), since 
it corresponds to differentiation in ^ under the change of coordinates £1 3 y — > 
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x(t, y) G Qlf, and we will work in both coordinate systems. This also avoids possible 
singularities in the change of coordinates. We will denote covariant differentiation 
in the ^-coordinates by V a , a = 0, ■ ■ • , n, and differentiation in the Xj-coordinates 
by di, i = 1, ■■■ ,n. The covariant differentiation of a (0, r) tensor k(t,y) is the 
(0, r + 1) tensor given by 

r)h 

\7 h - u ^- a r _ r d h T d k ; (9 A) 

y a^ax---a r — q L aa 1 hj d---a r J- a a r hj a, 1 ---di K^-^J 

where the Christoffel symbols T^ b are given by 

_ g cd ( dg bd dg ad dg ab \ _ dy c d 2 x i 

1 ab - TT + — ^ZI n..na~.h- \ A - C> ) 



2 V dy a dy b dy d J dx 1 dy a dy b ' 

If w(t, x) is the (0, r) tensor expressed in the ^-coordinates, then the same tensor 
k{t, y) expressed in the ^/-coordinates is given by 

k av .. ar (t,y) = —---—w h ... ir (t,x), x = x(t,y), (2.6) 

and by the transformation properties for tensors, 

_ dx % dx 11 dx lr dw h ... ir , , 

a ai -a r ~ dy a dyai " " Q y a r fai \ ' > 

Covariant differentiation is constructed so the norms of tensors are invariant under 
changes of coordinates, 

g a lbl . . .g-rbr^ arhi ^ = tfdl . ..^W n ... lr W n ... ]r . (2.8) 

Furthermore, expressed in the ^/-coordinates, 

d dy a d 
1 dx 1 dx 1 dy a 

Since the curvature vanishes in the ^-coordinates, it must do so in the ^/-coordinates, 
and hence 

[V o ,V 6 ] = 0. (2.10) 

Let us introduce the notation k a ... b ... c = g bd k a ...d... c , and recall that covariant differen- 
tiation commutes with lowering and rising indices: g ce V ' a kb-e-d = V a g ce kb- e ---d- Let 
us also introduce a notation for the material derivative 

d_ 

di 



*"5 



y=const 

Then we have, from j|, Lemma 2.2], that 

dx^ dx^ r f dv dv \ 

D ^-a r = Wl --- Wr [D tWiv ., r + ^-w,.. lr + ■■■ + —w n .., j . (2.12) 

Now we recall a result concerning time derivatives of the change of coordinates 
and commutators between time derivatives and space derivatives (cf. |4|, Lemma 
2.1])- 

Lemma 2.1. Let x = ft{y) be the change of variables given by ( |2.I|) , and let g a b be 
the metric given by ( p.2|) . Let Vi = 6ijV^ = v l , and set 

doc 1 

u a {t,y)= Vi {t,x) — , u a =g ab u b , (2.13) 



Then 



CHENGCHUN HAO AND T. LUO 

h ab =~A0a6, h ab =g ac h cd g db . (2.14) 

P dx l dx k dv l p dy a dy a dv k ^ 
<9?/ a <9y a &c fc ' dx l dx k dx % ' 

A<7a6 = V a u b + V b u a , D t g ab = -2h ab , D t dfx g = g ab h ab dfx g , (2.16) 

D t Y c ab = V a V b u c , (2.17) 

where d[i g is the Riemannian volume element on Q in the metric g. 

Proof. The proof is the same as that of 0, Lemma 2.1] except that we need to make 
some modification due to the difference of the definition of h ab . Indeed, the proof of 
( 2.15 ). ( 2.17 ) and the first part of (|2.16| ) is the same as the mentioned. The second 



part of (F1ED follows from (fTTJ) since = D t (g ad g dc ) = (D t g ad )g dc + g ad D t g dc 
(D t g ad )g dc + 2g ad h dc and then D t g ab = (D t g ad )5 b = (D t g ad )g dc g cb = -2g ad h dc g cb = 
—2h ab . The last part of ( f2.16j ) follows since in local coordinates d[i g = \f det gdy and 
D t {detg) = (det g)g ab D t g ab . ' □ 

We now recall the estimates of commutators between the material derivative D t 
and space derivatives and covariant derivatives V a . 

Lemma 2.2 ( [^j, Lemma 2.3]). Let di be given by ( |2.9| ). Then 

[D t ,d l ] = -{d l v k )d k . (2.18) 

Furthermore, 

[A,91=g-(^ 1 )(^)-9-, (2.19) 
where the symmetric dot product is defined to be in components 

and ^ r denotes the collection of all permutations of {1, 2, ■ • • , r}. 
Lemma 2.3 (cf. @, Lemma 2.4]). Let T ai ... ar be a (0,r) tensor. We have 

[D t , V a ]T ai ... ar = -(V ai V a u d )T da2 ... ar (V ar V a u d )T ai ... ar _ ld . (2.21) 

If A = g cd VcV d and q is a function, we have 

[D t , g ab V a }T b = - 2h ab V a T b - (Aw e )T e , (2.22) 
[D t ,V]q=0, (2.23) 
[D t , A]q = - 2h ab V a V b q - (Au e )V e q. (2.24) 

Furthermore, 

[A, V]g = E-( s + i) ( VS+1 ^) ■ ( 2 - 25 ) 
where the symmetric dot product is defined to be in components 

((V^r*) • V-q) ai ,„ ar = i £ ^Z^ a+ u d ) V d Z s+2 ... a „ r *- (2-26) 
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Proof. The proof is similar to that of |4], Lemma 2.4]. We only need to verify (|2.22|) 
and (|2.24|) since they involve the term D t g ab . Now from ( |2.16|) and (|2.21|) , it follows 
that 



[D t ,g ab V a ]T b =D t (g ab V a T b ) - g ab V a D t T b 

= (D t g ab )V a T b + g ab D t V a T b - g ab V a D t T b 
= -2h ab V a T b + g ab [D t ,V a ]T b 
= - 2h ab V a T b - g ab V b V a u e T e 
= - 2h ab V a T b - (Au e )T e . 
From (|1D and we have 



D t V a q =D t 

dx 



dx l 
dy a 



dig 



dx 1 
dy a 



dy a 
dx 



D t diq + di 
dv 

[D t , di)q + diD t q + dtq— 



dv 1 
dx i 



dx 1 



- (-diV k d k q + diD t q + diV £ d £ q) = —diD t q = V a D t q, 



Qya \ 1 I Qya 

namely, ( p^3D follows. Then, ( ggg ) follows from ( ggg ) and 

[A, A]g =D t Aq - AD t q = D t (g ab V a V b q) - g ab V a V b D t q 
= [D t , g ab V a ]V b q + g ab V a [D t , V b ]q 
= [D t ,g ab V a ]V b q. 
Therefore, we complete the proof. 
Denote 

dxi 



B t = SyB* = B\ a 
It follows, from ( |2.8|) , that 

\B\ 2 = \B\ 



B, 



(3 a = g ab (3 b , and \^ = f3 a /3 a . 



Qya- 



B, 



Qya 



(3a, P = P + 



8tt 



□ 



(2.27) 



(2.28) 



Then P = ^-s 2 on the boundary dQ. 



From flOp , (|T5a| ), ( ggg ), ( [TToD , Q, we have 



D t u a =D t vj 



dx j 



dxi 



dy a J dy 



-D tVj + VjD, 



dxi 
dy a 

dx k dvi 



dy 
-V Q P 



dx k 



1 dxi Q y b dy d dy c 

: d —PbO k — — V d Pc 



An dy a dx i dx 1 dxi 



ax J dx 1 



- V a P + — g bd g a eg ec PiN d f3 c + g oc u b V a u c 

47T 



47T 



/3 d V d /3 a + M c V a w c . 
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Similarly, we get 

„ dx j ^ _ „ dx j dx j „ dx k dv j 

D t Pa =-^—D t B i + BjD t — — = —B k d k v j + Ba— — ——- 
dy a 3 3 dy a dy a 3 3 dy a dx k 

=l3 d V d u a + (3 c V a u c . 

Thus, the system ( |1.1| ) can be written in the Lagrangian coordinates as 

1 



D t u a + V a P = u c W a u c + — /TV d /3 a , (2.29a) 

47T 



r = (/ \ ,,a r i — >" x 

D t p a = (3 d V d u a + (3 c V a u c , (2.29b) 

V a u a = in [0, T] x fi; V a /3 a = in {t = 0} x fi, (2.29c) 

\P\ = ? and /3 a iV a = on {t = 0} x (2.29d) 

p = on[0,T]xffi. (2.29e) 

3. The Energy Conservation and Some Conserved Quantities 

Firstly, the divergence free property of (3, i.e., div/3 = 0, is preserved for all 
times under the Lagrangian coordinates or in view of the material derivative, i.e., 
-Dfdiv (3 = 0. Indeed, from ( |2.22| ) and Lemma the divergence of (|2.29b ) gives 



D t (g ab V b (3 a ) = {D t ,g ab V b ]f3 a + g ab V b D t f3 a 
= - 2h ab V b (3 a - (Au e )(3 e + g ab V b ((3 d V d u a + (3 c V a u c ) 
= - 2h ab V b f3 a - (Au e )(3 e + V b (3 d V d u b + (3 d V d V b u b 

+ g ab \7 b (3 c \7 a u c + f3 c Au c 
= - g ac (V c u d + V d u c )g db V b f3 a + V b (3 d V d u b + g ab V b /3 c V a u c 
=0. 

Secondly, we assume that 

\Vu(t,y)\^C on [0,T]x9O, (3.1) 

then that (3 ■ N = is preserved for all times t in the lifespan [0, T], that is, we have 
13 ■ N = on [0, T] x dfl if (3 ■ N = on {t = 0} x dQ. Indeed, we have, from (|2~29"b"p 
and Lemmas |2.1| and [A.4| , that 

D t {(3 a N a ) =D t (g ab f3 a N b ) = N a D t (3 a + (3 a {D t g ab )N b + (3 a g ab D t N b 

=N a ((3 d V d u a + (3 d W a u d ) - V c u b f3 c N b - N d V d u a (3 a + (3 a g ab h NN N b 

=h NN /3 a N a , 

which implies, by the Gronwall inequality and the identity l-D^/H = \D t f\, that 

\((3 a N a )(t,y)\ <: e ct \(f3 a N a (0,y)\ = 0. (3.2) 

Thus, in view of the above three preserved quantities, the system ( j2.29j ), or fll.lQ , 
can be written in the Lagrangian coordinates as 

D t u a + V a P = u c V a u c + -^(3 d V d (3 a , (3.3a) 

D t Pa = P d V d u a + (3 c V a u c , (3.3b) 
V a w a = 0, V a (3 a = 0, in [0, T] x fi, (3.3c) 
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P = -^ 2 , \0\ = q, (3 ■ N = 0, on [0, T] x dQ. (3.3d) 
Finally, the energy defined by 

E (t) = J^\\u\ 2 + ±\^y^ g (3.4) 

is conserved. In fact, by (|2.16|) , ( |2.29|) , Gauss' formula and the fact D t d/i g = due 
to divu = 0, it yields 

jE {t) = J^D t (^-g ab u a u b + gV 6 AA) dfi g 

u a D t u a + i-/3 a A/3a) dfi g 
+ J ^(D t g ab ) (uaUb + ^PaPb) d^ 
[ [-u a W a P + u a u c V a u c + ^u a f3 d V d f3 a ]dfi g 
+ J (^(3 a l3 d V d u a + ^(3 a l3 c V a u)j dfi g 
- J h ab [u a u b + ^PaPbJ dfig 

N a u a Pdfi 1 + I u a u c V a u c dfi g + — I N d f3 d u a l3 a dH 
an Jn 4?r Jan 

+ 47 jjy^aUcdfig 

-\J 9 ac iScU d + V d u c )g db (u a u b + -^PaPb^j dfig 
=0. 

4. The First Order Energy Estimates 
From ( |2.21| ) and ( |3.3a| ), we have 
A(VftWa) + V 6 V a P 
= [A, V 6 ]u« + V 6 A«« + V 6 V a P 

= - (V a V 6M d K + ^V b ((3 d V d f3 a ) + V 6 KV a « c ) 

47T 

= - (V a V b u d )u d + — (V b f3 d V d (3 a + (3 d V b V d f3 a ) + V b u c V a u c + u c V b V a u c 

=V^ c V a w c + -^(V b f3 d V d (3 a + (3 d V b V d (3 a ). 
From and ( |3T3bD , we get 

A(V b /3 a ) = [A,V 6 ]/3 a + V b A/?a 
= - (V a V 6 n d )/3 d + V b (/3 d V dMa + /3 c V a u c ) 

= - (V a V 6 u d )/3 d + V 6 /? d V d u a + /3 d V 6 V d u a + V b (3 c V a u c + (3 c V b V a u c 
=V b (3 c {V c u a + V a u c ) + p d V d V b u a . 
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Thus, we obtain 

D t (V b u a )+V b V a P = V b u c V a u c + ^{V b f3 d V d (3 a + (3 d V b V d f3 a ), (4.1) 

47T 

D t (VM =V,/3 c (V c n a + V a u c ) + (3 d V d V b u a . (4.2) 

Now, we calculate the material derivative of g bd ^ ae V ' a u b V ' e u d . From Q2.16Q , (|2.14j ), 
(IQ5D , we get 

= (A^ M )7 ae V aMfe V e n, + 5 M (A7 ae )V«^V e! i d + 2g bd 1 ae (D t V a u b )V e u d 
= - 2g bc h cf gf d 1 ae V a u b V e u d - 2g bd 1 ac h cfl ^V a u b V e u d 

- 2g bd 1 ae V e u d V a V b P + 2g bd 1 ae V e u d V a u c V b u c 

+ ^9 bd l ae V e u d (V a (3 d V d f3 b + (3 d V a V d (3 b ) 

= - 7 ae (V c u/ + V } u c )V a u c V e u s - 2 7 a V e (V c M/ + V f u c )V a u d V e u d 

- 2 1 ae W e u b W a W b P + 2 1 ae W e u b W a u c W b u c 

+ ^^V e u b {y a p d V d f3 b + (3 d V a V d (3 b ) 
= - 2 1 ae V c u f V a u c V e u f - 4j ae ^ c V e u f V a u d V c u d + 2^ e V e u b V a u c V b u c 

- 2 1 ae V e u b V a V b P + ^i ae V e u b (V a p d V d p b + P d V a V d p b ) 

271 

= - 4 7 ae 7 /c V e M / V a M d V c M (i - 2 7 ae V e u 6 V a V 6 P 
+ ^-j ae VeU b (Va(3 d V d p b + f3 d V a V d f3 b ). (4.3) 

271 

Similarly, from (|4.2|) , we have 
D t {g bd j ae V a (3 b V e (3 d ) 
= - 7 ae (V cM/ + V f u c )V a (3 c V e P f - 2 7 a V e (V cM/ + V f u c )V a (3 d V e p d 

+ 2 1 ae V e f3 b (y a f3 c V c u b + V a /3 c V bMc + (3 d V d V a u b ) 
= - 2 7 ae V cM/ V a /3 c V e /3 / - 4 1 ac 1 fe V c u f V a f3 d V e f3 d 

+ 2 7 ae V e /3 b V a /3 c V cMfe + 2 1 ae V b u c V e f3 b V a (5 c + 2 7 ae (3 d V ' e /3 fc V d V \u b 
= _ Ar e l fc V e u f V a f3 d V c f3 d + 2 7 ae V e /3 6 V a /3 c V c « 6 

+ 2r e P d V a p b V d V e u b . (4.4) 
Thus, by combining ( |4.3|) with ( |4.4|) , we obtain 

A (V V e V a u b V e u d + ^g bd j ae V a f3 b V e (3 d 
= - 4 7 ae 7 fc V e u f V a u d V c u d - -i ae ^ c V e u f V a (3 d V c (3 d 

7T 

- 2V b (^ ae V e u b V a P - ^i ae fl b V e u d V a (3 d 
+ 2(V fe7 ae ) heU b V a P - ^l3 b V e u d V a (3 c 



--2(D t g ac )g bd (cur\u) ab (cur\u) cd + 4g ac g bd (D t V a u b )(cmlu) 
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+ -L 7 -V e n fe V a /3 d V rf /3 6 + ^i ae V e /3 b V a f3 c V c u b . (4.5) 

Now, we calculate the material derivatives of |curlw| 2 and | curl /3 1 2 . We have 
A|curlu| 2 =D t (g ac g bd {cuYlu) ab {curlu) cd ) 

i a -Cic „bdf n n - \ f i - 

)cd 

= - 2g ae gf c g bd {V e u f + V /Me )(curl M ) ab (curl M ) cd 
+ 4g ac g bd (cur\u) cd V a u e V b u e 
-Ag ac g bd (cm\u) cd V a V b P 

+ -g ac g bd (curlu) cd (Vaf3 e V e (3 b + (3 e V a V e (3 b ) 

71 

= - 4g ae g bd V e u c (curl u) ab (cur\ u) cd 

+ -g ac (cm\u) cd (V a (3 e V e (3 d + (3 e V a V e (3 d ). 

71 

Similarly, 

A|curl/3| 2 =2(D t g ac )g bd (cm\f3) ab (cm\f3) cd + 4g ac g bd (D t V a f3 b )(cm\ (3) cd 
= - % a V d V eM c (curl/3) af) (curl/3) cd 

+ Ag ac g bd ( C m\p) cd (V a p e (V e u b + V b u e ) + (3 e V e V a u b ). 

Thus, we can get 

A(|curlu| 2 + icurl/3| 2 ) 

= - Ag ae g bd V e u c (cmlu) ab ( C mlu) cd + -g ac (cm\u) cd V a (3 e V e f3 d 

7T 

- -^ M V eM c (curl/3) afe (curl/3) cd 

7T 

+ -g ac g bd (cur\p) cd V a (3 e (V e u b + V b u e ) 

71 

+ -V e {g ac (cm\u) cd f3 e V a f3 d ). (4.6) 

71 

Define the first order energy as 

Exit) =J ^g b V e V a u b V e u d + ^g bd r e V a f3 b V e {3?j dfi g 

+ J ^|curlu| 2 + ^-|curl/3| 2 ^ d\i g . (4.7) 

Let us recall the Gauss formula for Q and dQ: 

/ V a w a dfU g = / N a w a d f u J , and / V a f a d^ = (4.8) 
Jn Jan Jan 

if / is tangential to dQ and N is the unit conormal to dfl. 
Then, we get the following estimates. 



Theorem 4.1. For any smooth solution of MHD ( |3.3|) for ^ t ^ T satisfying 

|VP|<M, |V«| <M, zn[0,T]xfi, (4.9) 
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|0| + | Vu| + - ^K, on [0, T] x dQ. (4.10) 

We have for t G [0, T] 

E 1 (t) <: 2e CMt E l (0) + CK 2 (Vol SI + E (0)) (e CMt - l) . (4.11) 
Proof. By Q4.5Q, (|4.6|) and Gauss' formula, we have 



jE^t) = J^D t (g bd r e V a u b V e u d + -Lg bd ^V a (3 b V e (3 d ^j dfx g 
+ J D t |jcurlw| 2 + -^-|curl/3| 2 ^ dfi g 
+ f (g bd i ae V a u b V e u d + ^g bd ^ e V a (3 b V e f3 d \ihdfi g 



Jn \ 

+ J ^|curlu| 2 + -^|curl/3| 2 ) tr hd[i g 



-4 f 7 a V c V e n / V aM d V c ^ci/i 9 -- / j ae j fc V eUfVa^V^adfig 
Jn 71 Jn 

~ 2 J dn N b (l ae V e u b V a Pdn g - ^ ae f3 b V e u d V a {3?j rf/i 7 (4.12) 

+ 2 Jjybl ae ) (v e u b V a P - -^(3 b V e u d V a f3 d ^j dfig (4.13) 

+ 7^ / 7 ae V e u b V a (3 d V d p b dfi g + -L [ j ae V e p b V a p c V c u b d^ g 
^ Jn 271 Jn 

-4 / g ae g bd V e u c (curl u) ab (curl w) cd d/i 9 
+ - / £ ac (curl M ) cd V a /3 e V e /3^/i 9 



7T 



^ ae ^ M V eM c (curl/3) af) (curl/3) C(i (i / i s 

n 



+ Z / 9 ac 9 bd (cni\(3) cd V a p e (V e u b + V b u e )dfi g 
n 



71 



+ - f N e (3 e g ac (cm\u) cd V a /3 d d^. (4.14) 
n J on 



I on 

+ (g bd j m V a u b V e u d + -^g bd j ae V a (3 b V e (3 d j trhd^ 

+ / (|curlw| 2 H |curl/3| 2 J tr /id/i„. 

in V 47r / 

Since P = ^ 2 on <9fi, it follows that VP = 0, i.e., 7^V d P = 0, and then 
7 ae V a P = (? ce 7c V a P = on the boundary <9f2. In addition, f3 ■ N = on <9f2. Thus, 
the integrals in Q4.12Q and (|4.14j) vanish. 

From (|A.5|) and ( |A.3|) , we get 

^ = (^ - N a N c )V c N b = V a N b - N a V N N b = V a N b , (4.15) 
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since in geodesic coordinates VnN = 0. It follows that 

V 67 ae =V h {g ae - N a N e ) = -V b (N a N e ) = -(V b N a )N e - (V b N e )N a 



e^N e - e e b N a . 



Thus, by the Holder inequality, ( [4.1 0|) and Lemma |A.5| , we get 
lOI <CK (||V«|| i2(n) ||VP|| ioo(n) (Volfi) 1 ^ 

+ llVull^ m\ LHn) ||V/3|| i2(n) 

^CKM ((Volft) 1/2 + S 1/2 (0)) E\ /2 (t). 

For other terms, we can use the Holder inequality directly. It yields 
d 



j t Ei(t) ^CKM ((Volft) 1/2 + £ 1/2 (0)j E\ /2 (t) 



+ C||Vd 



L°°(n) 



(l|V^|| 2 L2(n) + ||V/3|| 



2 

L 2 (n) 



+ ll c urlw|| L2(n) + ||curl/3|| 

L 2 (Q) 



<:CKM ((Volfi) 1/2 + El /2 {0)j E\ ,2 {t) + CME 1 (t). 

From the Gronwall inequality, it follows that 

E\ /2 (t) ^ e CMt / 2 El /2 (0) + CK ((Volfi) 1 / 2 + £ 1/2 (0)) (e CMt ' 2 - l) , 

which implies the desired result. □ 

Remark 4.2. Since ( 4.12|) , especially the integral involving P, vanishes, we do not 
need the boundary integral in the first order energy Ei(t). But in higher order 
energies estimates, we need to introduce boundary integrals for P in order to absorb 
the analogy integral to ( }4.12|) . 

5. The General t-th Order Energy Estimates 
From d^TJ), (|2T9l ), (|L8a|), we get 



D t V r u a =D t V ai ■ ■ ■ V ar u a = D t 



dx 11 dx %r dx 1 



dx 11 


dx ir 


dx 1 


dy ai 


Qya r 


dy a 


dx 11 


dx lr 


dx 1 


dy ai 


Qya r 


dy a 


dx 11 


dx lr 


dx 1 


dy ai 


QyO, r 


Qyd 



Qyai Qya r QyC 



■ ■ ■ diVi 



D t d h ■ ■ ■ d ir Vi + jr—de ■ ■ ■ d ir Vi + ■ ■ 

OX n 

OX %T ox 1 



[D t , d r ]v t + d r D t Vi + p-d e ■ ■ ■ d tr v t + 

ox 11 



r-1 

E 



r 

s + l 



(d 1+s v) ■ d r - s Vi - d r diP 



1 

4tt' 



+ ^(B^M + ^-—d e ■ ■ ■ d lrVi + 



dx h 
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r-l 



! - s u a 



s=l ^ ' 

+ V a u c V r u c + ±±Q V s /3 c V r " s V c A 



=o 

where 

(v s /3 c v r - s v c /3 a ) ai ... ar = pvi q ..,jx;: i ....vA. (5-1) 

Thus, due to div (3 = 0, we get for r ^ 2 
AV r « a + V r V a P 

r-2 / x 

=(curl W ) ac W + sgn(2 -r^l/ ) (V 1+S n) • V r ' s u a 

s=l ^ ' 

+ ^V c (/3 c V^ a ) + ^ E Q V*/TV-V c /3 a , (5.2) 

s=l 

where sgn(s) is the signum function of the real number s, i.e., sgn(s) = 1 for s > 0, 
sgn(s) = for s = 0, and sgn(s) = —1 for s < 0. Of course, we use this notation 
sgn(2 — r) to indicate that the related term vanishes for r = 2. 
Similarly, by noticing that div /3 = 0, we have 

D t V r l3 a =V a u c V r f3 c - V r w c V c /3 a 

r-2 , v 

- sgn(2 -r)J2{ s r +1 ) (V 1+S u) ■ V r ~ s p a 

s=l ^ ' 
r 

+ V c ((3 c V r u a ) + E V s /3 c V r " s V cMa . (5.3) 

s=l 

Define the r-th order energy for r ^ 2 as 

^r(*)= / g bd l aS l AF V^ l V a u b V r F - l V f u d d^ g + I IV-Wlul 2 ^ 

+ 7^ / |V r - 1 curl/3|>,+ / -ft 1 AF V r ^V a PVf 1 V f P<&dn 1 , 
where $ = l/(— VjyP) as before. 

Theorem 5.1. Le£ r e {2, • • • , n + 1}, then there exists a T > suc/j £/iai the 
following holds: For any smooth solution of MED ( |3.3| ) for ^ £ ^ T satisfying 

\p\ <Mi /orr = 2, m [0,T] x fi, (5.4) 

|VP| ^ M, |Vu| <M, |V/3| < M, m [0,T] x fi, (5.5) 

|0| + 1/lq ^K, on [0,T] x dn, (5.6) 

-V^vP ^ e >0, on [0, T] x dQ, (5.7) 
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\V 2 P\ + \V N D t P\ ^L, 
we have, for t G [0, T], 

E r (t)^e Clt E r (0) + C 2 (e c ^-l) 



on[0,T]x<9Q, (5.8) 



(5.9) 



where C\ and G% depend on K , K\, M , M\, L, \je, Volfi, E (0), Ei(0), ■ ■ ■ , and 
£ r _i(0). 



Proof. 



Jt E ^ 



! D, (g^'l^V^W^Vf^ f u d ) in, 

JQ 

+ f D t (jVWV^V/fl,) dfig 

+ / D t \V r - 1 cur\u\ 2 dfi g + [ A|V r_1 curl/3| 2 d/i s 

+ / g^i^V^VaUbV^VfUdtrhdfig 
Jn 

+ / |V r_1 curlu| 2 tr/id// g + — / | V^curl /3| 2 tr 

+ / D t (i*f 1 AF V r A X V a FV r i 1 V f P)M lh 
Jan 

+ J^l^V^VaPV^VfP n±+trh- h NN ^j M^. 

We first estimate (|5.10|) , (|5.11 ) and ( |5.16 ) . From Lemmas [2.1| and [A.4| , and ( j5~2] 

we have 

=(A^ d )7 a/ 7 AF V^ 1 VaM 6 V r F - 1 V / n (i + rg^Dt-ff^ST^VaU^VfUi 

+ 2( ? b V / 7 AF A(V^ 1 V aMb )V^ 1 V /Md 
= - 2V c M e7 ^7 AF V^ 1 V aM c V r F - 1 V /M e - 4rV c « e7 a VV F V; 4 - 1 V f y V^V^, 

- 2 7 ^ 7 AF V^ 1 V /M f 'V^ 1 V a V fe P + 2 7 a V F V^ 1 V /M b (curl M ) 6c V^- 1 V aM c 



(5.10) 
(5.11) 
(5.12) 
(5.13) 
(5.14) 
(5.15) 
(5.16) 
(5.17) 



+ 2sgn(2 - r) 7 a/ 7 AF V r F - 1 V / ^ £ ( a + i) ((V s+ M ■ V r ~V) 



An 



Similarly, 

D t {g hd i af i AF V r A l V a p b V r f x V f p d ) 
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+ 2 1 a f 1 AF V r F - 1 V f (3 b V b UcV^ l V a l3 c - 2 7 a V F V^ 1 V / /3 fe V c /3 b V^ 1 V 

+ 2sgn(2 - r)7 a V F V r F " 1 V / /3 fc E ( s + i J (( vl+ ^) ' V r -ft) ^ 
+ 2 1 af l AF V J F - 1 V f (3 d V c ((3 c V r Aa u d ) 

r 

+ 2 1 af l AF V F - 1 V f (3 d J2 Q (V s /3 c V r - s V cMrf ) Aa , 



8=1 

and 



Thus, we get 

O + (^TTI) + O 

<C (l|VM|| LO c (n) + ||V/3|| LOO(n) )£; P (t) 



r-2 



+ CPy 2 (t) £ || V* 1 *^ (||V-«|| i4(n) + || V-'/?|| L4(n) ) (5.18) 

s=l 
r-1 

+ ^ /2 W E H^llx^n) (||V r - s+1 HL 4( n) + l|V-*-VL (n) ) (5-19) 

+ 2 7 a W^P (AV^P - JjV 6 V^j (5.20) 

+ 2 / V 6 ( 7 a/ 7 AF ) V'f^fvfVFVaPdfig (5.21) 

+ iV^V^-V/^V^/i, (5.22) 
2^ Jsn 

" 27 Su Vc ^VfUdFVdPdfig. (5.23) 

Due to /3 ■ iV = on <9£7, ( |5.22| ) vanishes. Let a be a (0, r) tensor and n G {2,3}. 
Then from Lemma |A.12| , we have, for i\ ^ l/Ki, that 

11/311^(0) <CJ2 < 2 " S HV s /3|| i2(n) ^(^)^ 8 1/2 (i)- (5-24) 

Thus, for the last integral, by the Holder inequality and the assumption (|5.5| ), we 
have for any r ^ 3 



O ^CK\\/3\\ Loom E r (t) ^ C{K,K X ) &E\f\t)\ E r 



(t). (5.25) 



For r = 2, we have to assume the a priori bound \/3\ ^ M\ on [0, T] x f2, i.e., (|5.4j) , in 
order to get a bound that is linear in the highest-order derivative or energy. Then, 
we have by ( |5.4|) 

( ggj ) ^ORT ||/3|| £oo(n) ^r(t) < C{K, M x )E r {t), for r = 2. (5.26) 
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By the Holder inequality, we have 

(H3)^L7^y 2 (t)||V r P|| L2(c) . (5.27) 

From ( 1.8a|) , we have 

dj(D t v j ) + AP = ^-d j (B k d k B*), 

47T 

which yields from fl2.18|) 

1 

Since A is invariant, we have 



AP = -d jV k d k v j + —djB k d k B j . 



AP = -V a u b V b u a + ^-V a f3 b V b (3 a . (5.25 



It follows that for r > 2 



V r - 2 AP =V r - 2 ( -V a u b V b u a + -^V a /3 fe V b /3 a 

47T 



" ( T 2 ) V s V a M fe V^ 2 - s V b u a 
s=0 V s / 

„_n \ / 



s=0 

From Q5.24j ), we have for s ^ 

2 2 

nv s /?iu (n) ^E< 2 ~1 v ^IIl^) ( 5 - 29 ) 

and, similarly, 

2 

HV s u|U (n) ^ C(^) J>$(*)- (5.30) 

1=0 

By Holder's inequality, ( |5.29| ) and ( |5.30| ), we get for r e {3,4}, 

||yr-2 a nil 

r-2 r-2 

^ E 1 1 V s V a u b V r - 2 - s V b u a | L, (0) 1 1 V s V a (3 b V r - 2 ~ s V b p a 1 1 L2(n) 

T" — 1 „ . 1 1 i /"HIV7/QII llwr— 1, 



^ ||V«|| ioo(f2) || V-^H^^ + C7 || V0|| LOO(n) || v- 1 /? 

+ (r - 3)C7 (||V 2 M || LOO(n) ||V 2 n|| L2(n) + || V 2 /3|| ioo(n) ||V 2 /3|| 



R\\ 

L 2 {VL) 



^C(K x )E¥&)Y^E\l*{t) + (r - 3)C7(P 1 )^2 1/2 (^) E^ /2 ^) 

£=1 fc2 

r-1 

^C(tfi) Et(t) + L7(P!)P 2 1/2 (t)Py 2 (t). (5.31) 



1=1 
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, i2(n) ||V«|| L2(n) ||V«|| LOO(n) + C || V(3\\ L2(n) ||V/3|| Loo(n) 



For r = 2, we have a simple estimate from the assumption (|5.5|) and Holder's 
inequality, i.e., 

I API 

^CME\^{t), (5.32) 

which is a lower energy term. Thus, by flA.,17 ), Q5.31 ) and ( [5.32D , we obtain for any 
5 r > 

l|V r P|| i2(f2) ^S r \\UV r P\\ L2(m) 



+ C(l/5 r: K,Vo\n) II V s AP| 



L 2 (Q) 



s<r-2 



r-1 



nnv r p 



L W + CWr, K, K h M, Vol Q) ^ P,(t) 

£=1 



+ (r -2)C(l/5 r ,K,K 1 ,M,Vo\tt)El /2 (t)El /2 (t). (5.33) 

Now we estimate the boundary terms. Since P = ^ 2 on <9f2, by ( |A.18|) , we have 
for r > 1 



nv r P|| LW ^(^^(ll^l 



L°°(9f2) 



2 > £ IK 



L 2 (dQ) 



X 



E ll v " p l 



L 2 (9Q) ' 



(5.34) 



fe<T--l 



From O, we get nV 2 P = tfVjyP and then, by Q, ( FBI) , ( CTD , Q and 
(|5.33| ), we get 



\L 2 (an) 



nv 2 p 



v^p 



L 2 (an) 



^ - ||nv 2 p| 

e 11 1 



L 2 (an) 



(5.35) 



nv 2 p 



L 2 (dU) 



\L°°(an) 11^* Hz, 2 (an) 



^C(^,Volfi)(||V 2 P|| L2(n) + ||VP| 



L 2 (n) 



^c{K,Vo\n)5 2 ||nv 2 p| 



L 2 {dn) 



C(K,Vo\Q)(yo\Q) 1/2 M 



+ C{1/5 2 ,K,K 1 ,M, Vol fi)Pi(t), 



(5.36) 



where the first term of the right hand side of ( |5.36| ) can be absorbed by the left 
hand side if we take 5 2 so small that, e.g., C(K, Vol Q)5 2 ^ 1/2. Thus, it follows 
that 



|nv 2 p| 

72 



\ L2(m) ^C(K,K 1 ,M,V61Q)(1 + E 1 (t)), (5.37) 
! V2p L 2 (n) <C(K,K l7 M,Vo\n)(l + E^t)), (5.38) 
I^IIl 2 (^) ^C , (K,K 1 ,M,Volfi,l/£)(l + P 1 (t)). (5.39) 
By Theorem |4.1| , there exists a T > such that Ei(t) can be controlled by the 



initial energy £?i(0) for t E [0,T], e.g., E x (t) < 2Pi(0). Thus, from flOl , flOSl) , 
(|5.5| ) and (|5.38|) we have 



l nv3p IL 2 (^) &{K, K,) {K + ||0|| i2(9fi) ) £ II V*P| 

fc<2 



Z 2 (<9Q) 
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^C(K, K u M, Vol Q, l/e)(l + E^t)) ^ || V fc P| 



\L 2 (n) 



k^3 
73 ; 



<C(#,tfi,M,Volfi,l/£,.E7i(0)) ||V 3 P| 
+ C(AT, Ki, M, Vol n, l/e, Pi(0)). (5.40) 

From (|533|) , 

ll v3p L ( n) Wifi.M^olfi,!/^!^)) ||V 3 P|| i2(n) 
+ (5 3 c7(ir, Ki, M, Vol Q, l/e, Ei(0)) 

2 

+ C(l/<f 3 , # , M, Vol Q) Et(t) 

i=i 

+ C(l/5 3 , K, K u M, Volfi)P 2 1/2 (t)P 3 1/2 (t), (5.41) 
which, if we choose £3 > so small that 

5 3 C7(^,iif 1 ,Af J Voin,l/e,Ei(0)) < 1/2, 

yields 

2 

II v3p L(«) ^> M ' Volfi ' Ve, ^i(O)) + C(K, K u M, Volfi) ^ £*(t) 

£=1 

+ L7(iT, Kx, M, Vol ft)P 2 1/2 (t)P 3 1/2 (t), (5.42) 

and then 

ll nv3p IUo» < C ( K > K ^ M ' v ° m > Ve> ^1(0)) 



X 



l + ^^) + ^2 /2 (^3 1/2 W • (5-43) 



Since 

V b V N P = 7b d V,(iV a V a P) = (5t - AyV d )((V,iV a )V a P + V a V d V a P) 

=ey a p + v a v 6 v a p - N b N d (e a d v a p + v a v d v a p), 

from ( A.31|) , it follows that 
|| VVjvP|| i2(an) ||0|| L oo (a n) ll VP llL2(an) + C II v2p IIl2(M) 

<c(#, voin) ( ||v 3 p|| L2(n) + ||v 2 p|| i2(f7) + ||vp|| 

2 

<C(tf, fsfi, M, Vol ft, l/e, Pi(0)) + C(tf, Ki, M, Vol ft) ^ 

+ C(K, K u M, Vol ft)P 2 1/2 (£)P 3 1/2 (t). 
Thus, by flOD , it follows that {V6)V N P = nV 3 P - 3#®VVjvP and 

IFi i2(an) ^ (||nV 3 P|| L2(an) +C||0|| Loo(aj) ||vv^p|| l2(9Q) ) 
id, M, Volfi, l/e, Pi(0)) 

l + E^W + P 2 1/2 (t)P 3 1/2 (t)j. (5.44) 



22 CHENGCHUN HAO AND T. LUO 

Hence, from flOp , flOU) , it yields 

ll n ^L 2(9Q) <c(k, k x ) (k + \\e\y m + 1| ve\\ L2[m] ) 

Then, from ( |5.33| ), we can absorb the highest order term || V 4 P|| l2 (q) by the left 
hand side for 64 > small enough which is independent of the highest energy E^(t), 
and get 

II v4p L ( n) ** C ( K > K ^ M ' Vo1 ^ l l £ i p i(°)) 

x ^1 + f^E e (t) + E¥ 2 (t)El /2 (t)j , (5.46) 

ll nv4p L ( ^) < C ( K > K ^ M ' Vo1 ^ l / £ i E M) 

x ^1 + f^E e (t) + El /2 (t)El /2 (t)j . (5.47) 

Therefore, from (|5.38|) , (|5.42| ) and (|5.45|) , we obtain for r ^ 2 
\\V r P\\ L2{n) ^C(K, K x , M, Vol J], 1/e, E 1 (0)) 

' r-l \ 

1 + J>,(f) + (r - 2)El /2 (t)El/\t) , (5.48) 



which, from ( |5.27| ), implies 

(CT) <C(if, Ki, M, Vol 0, A(0))£ r 1/2 (t) 



r-l 



1 + 5>,(t) + (r - 2)P 2 1/2 (t)Py 2 (t) . (5.49) 
1=1 / 

Now, we turn to the estimates of ( J5.20 ). Since P = ^<j 2 on dQ implies 7^ V a P = 
on dQ, we get from ( |A.3| ), by noticing that $ = — 1/VVP, that 

jV 6 =V N PN b = N a V a PN b = 5y a P - 7b a V a P = V fe P. (5.50) 

By the Holder inequality and ( |5.50| ), we have 

(PI ^ n*ui2 (fln) p r 1/2 (t) ||n (a (vp) - tr^vV) || L2(an) 

=c Il^lli2 (fln) £ r 1/2 (*) ||n (a (v^p) + v r u ■ vp)\\ LHdn) . (5.51) 

By ( g25D , it follows that 

D t V r P+V r u ■ VP = [A, V r ]P + V AP + V'u ■ VP 



r-2 , s 

sgn(2 - r) ^ I r J (V s+1 m) • V r ~ s P + V AP 

s=l ^ ' 



(5.52) 



We first consider the estimates of the last term in (|5.52| ). By ( |A.18[) and ( |A.31| ), 

we get, for 2 < r < 4 



||nv r AP|| 



L 2 (dn) 
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^C(K, K u Vol tt) M|t/|| LtMW) 

xEll vfc ^ p W 



+ (r-2) || V ^ 



L 2 (dn) 



(5.53) 



From (| A . 1 7| ) , it follows that 
l|V r AP|| L2(C ) 

<6\\nv r D t p\\ L2{dn) + c(i/5,K,Vom) J2 HV s AAP|| i2(n) 



(5.54) 



s^r-2 



By (gjp , flOp , Lemma [D] (Q, Q and fl3j), it yields 
AAP =2/i a6 V a V fe P + (Au e )V e P - D t (g bd g ac V a u d V b u c ) 

=2/i ab V a V b P + (A« e )V e P - 2A(5 M )V B «dV 6 « fl 

- 2( ? M A(V a n,)V b n a + ^-D t (g bd )g ac V a p d V b p a 

+ -Lg bd D t (V a p d )V b (3 a 
=2h ab W a W b P + (Au e ) V e P + 4/i M V a u d V 6 u a - i^V a /3 d V b /3 a 

7T 

+ 2g bd V b u a V a V d P - 2g bd V b u a V a u c V d u c 

- ^ b u a (V a (3 c V c (3 b + l3 c V a V c (3 h ) 

+ ^9 bd V b f3 a (V a f3 e (V e u d + V d u e ) + f3 e V e V a u d ) 

Z7T 

=4g ac V c u b V a V b P + (Am 6 ) V e P + 2V e u b V b u a V a u e 

- ^v bU a V a (5 c V c (3 b - ^v bU a (3 c V a V c (3 b + -Lv b f3 a (3 e V e V a u b . 

Z7T Z7T Z7T 



By (15729D , (|5T33|) and Lemma [ATT21 , it follows that for s ^ 2 
||V s AAP|| i2(f ,) 

<C \\Vu\\ Loo{n) \\V s+2 P\\ L2{n) + s(s - l)C || V 3 u\\ L2{n) ||V 2 P|| LOO(Q) 



+ sC W^Lnn) II v s+1 P\\ LHn) + c || ^u\\ L2(n) \\VP\\ Loa{n) 



l|V/3|| Loom) ||V/3|| 



Il 4 (h) 



lL°°(n) 

IV 2 w| 



L°°(Q) 



lL 4 (Q) 



+ c(||v 

+ S ( S -l)L7||V M || Loo(n) ||V 2 M 
+ C\\Vu\\ Loa{n) ||V/3|| 

+ sC ll^lUn) II^Ljn) ((« - 1) l|V/3|U (n) + 
+ 8 (a - 1)C \\Vu\\ Lao{u) \\V 2 f3\\ LHn) \\V 2 (3\\ L4(n) 
+ C\\Vu\\ Lc 



7S+1 



L 2 {il) 



\L ll^ S+1 ^llL2(f7) 



L°°(f2) 



'(«) 



II V s+2 /3ll 
L°°(n)|| v ^llL2(n) 



1/2 
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+ sC \\ V Hmn) WWl-v) (( s - 1) 11^11^(0) + HV/3|U (n) ) 
+ s ( s - 1)C ||V 3 /3|| L2(n) ||/?IU (n) ||V 2 M || ioo(n) 
+ s ( s -l)CjV/% 0O(n) ||V 2 /3|| 

L4(Q) H^HIl^) 

+ s ( a - l)C \\VP\\ Loo(sl) WWi^w || V 4 «|| L2(n) 

+ s(s - l)C ||V 2 /3|| LOO(n) W\\ L ^ a) || V 3 «|| L2(n) . (5.55) 
From ([ATTD , flOO)) , it follows that 

l|V s+1 HL, ( ,) <C\\V°u\til m (i2\\V^u\\ L2(n) Kt 

\e=o 

2 

^)X;^S(t). (5-56) 

We can estimate all the terms with L 4 (fi) norms in the same way with the help of 
( 5.29j ). ( |5.80| ), the similar estimate of P and the assumptions. Thus, we obtain the 
bound which is linear about the highest-order derivative or the highest-order energy 
E l r /2 {t), i.e., 

|| V s AD t P\\ L2{n) ^C(K, K u M, Mi, L, 1/e, Vol ft, E {0)) 

r-l 

x (l + ^^W)(l + ^ 1/2 W)- (5-57) 

1=0 

Thus, from (|5.53 ), Q5.54Q , (|5.57 ) and taking some small 5's which are independent 
of E r (t), we obtain, by induction argument for r, that 

\\UV r D t P\\ L2(m) ^C(K, K h M, Mi, L, 1/e, Volfi, E (0)) 

r-l 

x (l + £i^)(l + £ r 1/2 (*)). (5.58) 

1=0 

To estimate Q5.52j ), it only remains to estimate 

||n ((V s+1 u) ■ V r ~ s P) \\ L2{m) for 1 ^ s ^ r - 2. (5.59) 

For r = 3, 4 and s = r — 2, we have, by ( |5.8| ) and Lemma [A.14j , that 

||n((v^)-v 2 p)|| L2(an) 

< IIV r-1 wll IIV 2 Pll 

^ II v u \\L 2 (dn) II v r \\L°°(dn) 



| L 2(n-l)/(n-2)( 9 Q) 

. -L- IIW- 



^CL \\V 2 u\ 

^C(K,Vo\Q)L (||V r M|| L2(f7) + ||V r-1 w| 
^C{K,L, VolO) (^(t) + P r 1/2 (t)) • (5.60) 
For n = 3, r = 4 and s = 1, by ( |A.6| ), Lemma [A.14| and ( [5.33| ), we get 



|n ((V 2 «) • V 3 P) 

|nv 2 w • nv 3 p + n(v 2 w ■ n)®ti(n ■ v 3 p) I 



L 2 {dn) 



L 2 (dn) 
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< c ll nv2 «L ( an) ll nv3p L 4( an) 



+ L7||n(^v 2 Ma )|| L4(m) ||n(v^v 2 p) llL4(9C) 



||V 2 n|| L4(QC) ||V 3 P|| i4(sc) 



^C(K, Volft) (||V 3 n|| L2(n) + ||V 2 M || L2(n) ) (||V 4 P|| L2(Q) + ||V 3 P| 



L 2 (n) 



3 



<c(#,#i,Voift)(£ 3 1/2 ^ + p 2 1/2 (t)) £>.(*) + $>i /2 (t) A 1/2 W 

\s=0 \J=0 / 

3 4 

^L7(iT, K 1; Vol ft) P s (t) ^ P f 1/2 (t). (5.61) 

s=0 1=0 

Hence, we have 

(IQOD ^C(ir, i^ x , M, Mi, L, 1/e, Vol ft, P (0)) 

r-1 

x (l + ^P s (t))(l + P r (t)). (5.62) 

s=0 



By Lemma A.ll, we can obtain 



r-1 



( gig) + flOp sC C(K,K h M, Vol ft, 1 + ^P,(t) P r (t). (5.63) 



s=0 



Therefore, we have shown that 

CT ) + ( pTTD + (|515D ^C(K, K 1; M, Mi, L, 1/e, Vol ft, P (0)) 



r-1 

We now calculate the material derivatives of | V r_1 curl w| 2 and | V r_1 curl/3| 2 . From 



i — i 

(l + ^P s (t))(l + P r (t)). (5.64) 



Lemma |2.1| , fl5.2|) and ( |5.3| ) , we have 
A ( |V r -Wd 2 + — |V r -W/3| : 

\ 47T 

=A {g ac g bd g AF V^\ C m\u) ab V r f\ C m\u) cd ) 
+ i-A (^V F V^ 1 (curl/3) ab V r F - 1 (curl/3) cd ) 

= (r + l)D t {g ac )g hd g AF V r A\cuT\u) ab V r f\cm\u) cd 
+ 4^W F A (Vr'Vaiij) V r f\ C m\u) cd 

+ ^A(^ c )rt AF V r A - 1 (curl/3) afc VJr 1 (curl/3) cd 
+ —g ac g bd g AF D t (V^ x V a /9 6 ) V^curl/S)*, 

7T 

= - 2(r + l)^ ae V e n c ( 7 & V F V^ 1 (curln) a;) V r F - 1 (curln) c , 



^^ e V e nW F V^ 1 (curl/3) a6 V^ 1 (curl/3) 

,ac fcrf AFwr— 1 



cci 



- Ag ac g bd g A1, W r F {c\xi\u) cd S7 r Aa W b P (this vanishes by symmetry) 
+ ^g w g AF V r f\cm\ n) c ,(curl n) 6e V^ e 
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+ 4sgn(2 - r)g ac g AF ^ F -\ C m\u) cd E ( s + i J (^ 1+s u) ■ V r ~ s u d ) Aa 

s=l 

+ -sgn(2 - r)g ac g AF V r f\ C m\p) cd Y, ( J ((V 1+ %) ■ V r " s /3 d ) Aa 

+ -g™g bd g AF V r f\cv<iipUV r A J e V b u e - -g ac g AF V r f\ C m\(3) cd W r Aa u e V e (3 d 

7T 7T 

+ -V e {g ac g AF W r F -\ C m\u) cd V A J d ) 

71 



1 - 

+ -^^^curlttj^J] Q (V s /3 e V r " s V e /3 d ) 

s=l 



Aa 



Aa ' 



Noticing that (3-N = on <9f2, then by the Holder inequality and the Gauss formula, 
we get 

(|512j) < C(K, K u M, Volfi, 1/g) (l + E^W) ^r(t). (5.65) 



s=0 



Thus, by ( |Q3 ) and (ggg) , we get 

A(VatP) =D t (N a V a P) = (D t N a )V a P + N a D t W a P 
= (-2h a d N d + /i^iV") V a P + N a W a D t P 
= - 2h a d N d V a P + h NN V N P + V N D t P, 

which yields 

& t D t V N P 2h a d N d V a P V N D t P 

J = "v^P" = hNN + "V^' (5 " 66) 

Thus, we can easily obtain that the remainder integrals, i.e., (|5.13|) , ( |5.14|) , (|5.15 ) 
and (15171) , can be controlled by C(K, M, L, l/e)E r (t). 
Therefore, we obtain 

j t E r (t) ^C(K, K l} M, Mi, L, 1/e, Volfi, P (0)) 

r-1 

X 

s=0 

which implies the desired result ( |5.9| ) by Gronwall's inequality and the induction 
argument for r £ {2, • • • , n + 1}. □ 



r— l 

(l + 3>(*)) (! + £»■(*)), (5-67) 



6. Justification of A Priori Assumptions 

Let JC(t) and e(t) be the maximum and minimum values, respectively, such that 
(|5.6|) and ( |5.7| ) hold at time t: 

/C(t)=max(||^,-)|| Loo(m) ,lAo(t)), (6.1) 
£(t) =||l/(V^P(t l -))||i i - (fl0 ) = l/e(t). (6-2) 
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Lemma 6.1. Let K\ ^ be as in Definition \A.$ , £(t) as in (|6.2|) . JTien there 
are continuous functions Gj, j = 1, 2, 3, 4, such that 



v u \\L™(n) 

+ 11^/311^(0) + 



li°°(0) 

|VP|| LOO( o) + ||V 2 P 



L°°(fi) 



^^2(^1, £,E Q ,--- , i£ n+ i, VolO), 
I ^ IIl<» (an) ^.Gz{K\,8,E Q , - • • , E n+1 , Volfi), 
VA^IL-ram <G 4 {K U £,E Q ,- ■■ , E n+1 , Volfi). 



(6.3) 
(6.4) 
(6.5) 
(6.6) 



l£°°(an) 

Proo/. follows from fl£5Cp , (|5~29| ) and flOl - From Lemmas and [O0| , 

we have 

2 

IIVPI 



lioo(n) ^(^oEII^^I 



L 2 (n) ' 



£=0 

n+l 



^IU {fln) ^C(*i)£||Vp| 



L 2 (9C) 



(6.7) 
(6.8) 



Thus, (0) follows from 0, (§3), Lemmas |A~T3T - TQ31 , ( gggj ), ( CT) and (^). 
Since, from (|A.7|) , 

|v 2 p| ^ |nv 2 p| = \v N p\\e\ ^ s-^ei (6.9) 

so (|6.5| ) follows from ( |6.4|) . ( |6.6| ) follows from Lemma |A.10| , ( |5.54|) , ( |5.57|) and 



Lemma 6.2. Lei ^ 1/ii and £1 6e as in Definition \A.3j . Then 

r 

^ C r {K x , £, Eq, ■ ■ ■ , E n+1 , Vol Q) ^ E s , (6.10) 



d F 



s=0 



and 



dt 



^C r (K u £,E ,--- ,P n+1 ,VolO). 



(6.1L 



Proof. ( |6.10| ) is a consequence of Lemma [TT] and the estimates in the proof of The- 
orems |4.1| and |5.1| . ( |6.11| ) follows from 

2 

\\V N D t P(t, -)llicx>( S n) 



d 


1 




^ c 


1 


dt 


-V N P(t,-) 




-VatP(V) 



L°°(en) 



and (KQ). 



□ 



As a result of Lemma |6.2| , we have the following: 

Lemma 6.3. There exists a continuous function T > depending on K\, £(0), 
E o (0), ■■■ , E n+1 (0), Volfi such that for 

00 ^rc^^o),^),-- - ^n+i(o),voin), (6.12) 

t/ie following statements hold: We have 

E 8 (t) ^ 2P S (0), + £(t)<2£(0). (6.13) 

Puri/iermore, 

^(o,y)y a F 6 



<^(t, 2 /)F a F b ^2^(0,|/)F a F b 



(6.14) 



28 CHENGCHUN HAO AND T. LUO 

and with e\ as in Definition \A.3j , 

\Af(x(t,y))-Af(x(0,y))\ <; 

\x(t,y) - x(t,y)\ < 



.£1 
16 : 

16 : 



9ac(t,y) 9(0,?/) 



dy 



dy 



16' 



y g an, 
y e ft, 
y g an. 



(6.15) 
(6.16) 

(6.17) 



Proof. We get fl6T3D from Lemma |6j if T(#i, £ (0), £ (0), • • • , £n+i(0), Volft) > 
is sufficiently small. Then from (|6.13|) and Lemma we have 



|V 2 P 



V«|| LO c (n) + l|V/3|| £0 o ( n) + l|/3|L« C n) + II vp IIl~(q) 
<C(lfi,£(0),So(0),... ,K+i(0)), 

+ 11^ 



< C{K u £{0), E (0), • • • , K+i(0), Vol ft), 
||VAP|| Loo( an) ^C(^i^(°)^o(0),--- ,K + i(0),Volft). 

By and (|OD , we have 

|AV«| < |V 2 P| + |Vu| 2 + |V/3| 2 + \/3\ | V 2 /3| , 
|AV/3| < |V/3| |Vu| + |/?| |V 2 w| . 
By ( |Q5D , flOID , Lemma |0| and (gig ), we have 

l|V«|| £ ac (an) + ||V0|| ioo(an) ^C(A^,£(0),P (0), ••• ,K+i(0),Volft). 



(6.18) 

(6.19) 
(6.20) 

(6.21) 
(6.22) 



Thus, by noticing that \/3\ = q on 5ft, it follows, from fl6.18| ), ( |6.19| ), Lemmas |A.1Q 
and |Q1 , QOOl ) and ( P9D , that 



|| AVM|| L0O(a0 ) 
< I|V 2 P" 



lL°°(8n) 



d V 2 M 



AV/3|| ioo(M) 

+ (l|Vw|| £ oo (an ) + l|V0|| L oo (an ) 
'|V 2 /3 



\L°°(dn) 



^C(K U £(0), P (0), • • • , P„+i(0), Vol ft) 



x 1 + || Vd 



l|V/3|| 



L°°(ar2) 



which yields, with the help of Gronwall's inequality, for ^ t ^ T 
l|V^,-)IU (8n) + ||V/3(t,-)|| L 
<e o(^(D)^(D),..^ l( a),vc l n)t (||Vt*(0, 011^ + ||V/3(0, .)L- (fln)/ 
_l_ e C(ftr 1 ,£(o),Eo(o),-,Bn + i(o),Vom)i _ ^ (6.23) 

If T is sufficiently small, it follows, after possibly making T > smaller, that 

||v«(r, Oilman) H-Iiv^r, OIL 

(6.24) 



which also guarantee the a priori assumption of ( |3 . 1| ) . 
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By (ET23p , QOD , (IQj ), (|5T57D and ( ^58|) , we have 
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IAVP|| LOO(n) = ||VA^IIl- (0 ) < C(^i) E l|v m AP| 



L 2 (Q) 



£=0 



^(i^, £(0), E (0), • • • , E n+1 (0), Vol 0), 
which implies for sufficiently small T > 

l|VP(t,-)IU(n)^2||VP(0,-)IU (n) . 
By and (gig ), we have 

< l|VP|U (n) + ||/?|| L « (n) ||V/3|| Loo(n) 

^(^^(O),^),-'- ,^n+l(0)), 



which yields 



K*,Oll i - (at) <2||T;(0,.)|| £oo(n) . 



(6.25) 

(6.26) 
(6.27) 
(6.28) 

(6.29) 



( |6.14[ ) follows from the same argument since D t g ab = V a u b + V b u a and by ( |6.18| ) 



\g ab (T,y)Y a Y b -g ab (O iy )Y a Y b \ < / | D t9ab (s, y)\ dsY a Y b 

Jo 



^2 



1 



a\rb 



\ Lac(n) dsY a Y b ^-g ab (0,y)Y a Y 



|V a M fe (s) 

if T is sufficiently small. Now the estimate for M follows from 

D t n a = h NN n a , 
and the estimates for x and dx/dy from 

D t x(t,y) =v(t,x(t,y)), 

dx dv(t,x(t,y)) dv(t,x)dx 



(6.30) 
(6.31) 



D t — = = 

dy dy dx dy' 

and (|6.29p and Q6.24 ), respectively. 

Now we use (|6.14|) - (|6.17|) to pick a Ki, i.e., l±, which depends only on its value 
at t = 0, 



(6.32) 
(6.33) 
□ 



Lx(t) ^ ii(0)/2. (6.34) 
Lemma 6.4. Let T be as in Lemma \6.% Pick i\ > sitcA i/iai 

1^(0, y x )) -N(x(0,y 2 ))\ ^ j, whenever \x(0, y x ) - x(0, y 2 )\ ^ 2i x . (6.35) 

Then ift^T, we have 

\J\f(x(t, yi)) - J\f(x(t, y 2 ))\ < £i, whenever \x(t, yx) - x(t, y 2 )\ ^ 2ii. (6.36) 

Proof. ( |6.36| ) follows from fl6.35f) , (|6.15|) and ( |6.16| ) in view of triangle inequalities. 

□ 



Lemma |T4| allows us to pick a K 1 depending only on initial conditions, while 
Lemma gives us T > that depends only on the initial conditions and K\ 



such that, by Lemma |6.4| , ^ K\ for t ^ T. Thus, we immediately obtain the 
following theorem. 
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Theorem 6.5. There exists a continuous function T > such that if 

T^T(/C(0),£(0),£ (0),--- ,£n+i(0),VolO), (6.37) 

any smooth solution of the free boundary problem for Euler's equations ( |1.1| ) and 
(PI ) forO^t^T satisfies 

n+l n.+ l 



E s (t) «C 2 # s (0), ^ t < T. (6.38) 



s=0 s=0 

Appendix A. Preliminaries and Some Estimates 

Let N a denote the unit normal to dtt, g a bN a N b = 1, g ab N a T b = if T G T(dfi), 
and let N a = g a bN b denote the unit conormal, g ab N a Nb = 1. The induced metric 
7 on the tangent space to the boundary T(dQ) extended to be on the orthogonal 
complement in T(Q) is then given by 

lab = g ab - N a N b , r" = 9 ab -N a N b . (A.l) 

The orthogonal projection of an (r, s) tensor S to the boundary is given by 

(ns)£X = ill ■ ■ ■ • ■ ■ <sz:: c i, (a.2) 

where 

j c a = 8% - N a N c . (A.3) 
Covariant differentiation on the boundary V is given by 

VS = nV5. (A.4) 
The second fundamental form of the boundary is given by 

e ah = (itvau = iy c N b . (a.5) 

Let us now recall some properties of the projection. Since g ab = 7 ab + N a N b , we 
have 

U(S ■ R) = U(S) ■ U(R) + U(S ■ N)®U(N ■ R), (A.6) 

where S®R denotes some partial symmetrization of the tensor product 5* <g> R, i.e., 
a sum over some subset of the permutations of the indices divided by the number of 
permutations in that subset. Similarly, we let S~R denote a partial symmetrization 
of the dot product S ■ R. Now we recall some identities: 

nV 2 g =V 2 q + 6V N q, (A.7) 

nV 3 g =V 3 g - 26®(erVq) + (V6)V N q + 36®W N q, (A.8) 

IIV 4 g =V 4 g - 6® ^5(W)-Vg + 8^V 2 g) - 2(V6)®(6~Vq) 

+ (V 2 9)V N q + 4{V6)®VV N q + 66®V 2 V N q 

- 39®(F6)V N q + W®9V 2 N q. (A.9) 

Definition A.l. Let M{x) be the outward unit normal to dS> t at x G d%. Let 
dist (xi,x 2 ) = \xi — x 2 \ denote the Euclidean distance in M. n , and for X\,x 2 G dS> t i 
let dist g@ t (xi, x 2 ) denote the geodesic distance on the boundary. 
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Definition A. 2. Let dist (x, d3>t) be the Euclidean distance from x to the boundary. 
Let to be the injectivity radius of the normal exponential map of dS>ti i- e -; the largest 
number such that the map 

d% x (-t , t ) ->• {x G R n : dist (x, d%) < t} 

given by (x, t) — > x = x + iN{x) ^ 

is an injection. 

Definition A. 3. Let < E\ < 2 be a fixed number, and let ti = ti(ei) the largest 
number such that 

\Af(x~i) — Af(x 2 )\ ^ £\ whenever \x\ — x 2 \ ^ t 1; Xi,x 2 G d@ t . (A. 11) 

Lemma A. 4 ( |4], Lemma 3.9]). Let N be the unit normal to dVt, and let h a b = 
\D t g a b- On [0,T] x dVL, we have 

D t N a = h NN N a , D t N c = -2h c d N d + h NN N c , (A.12) 

D tl ab = -2 7 ac /w 7 d6 , (A.13) 

where h^N — h a bN a N b . The volume element on dQ satisfies 

Dtdfi^ = (trh- h NN )d^ = (tr 9u ■ N + >y ab V a u b )d^, (A. 14) 

where Ub denotes the tangential component of Ub to the boundary dQ. 

Lemma A. 5 (cf. @, Lemma 5.5]). Let w a = w Aa = V r A f a , V r A = V ai - • • V ar , / be 
a (0, 1) tensor, and [V a , V&] = 0. Let divw = V ' a w a = V r div/ ; and let (curlw) ab = 
V a w b - V b w a = V r (curl/) afe . Then, 

\Vw\ 2 ^ C(g ab Y d l AB V c w Aa V d w Bb + |divw| 2 + |curluf ). (A.15) 

Lemma A. 6 ( [§], Proposition 5.8]). Let t and i\ be as in Definitions \A.t\ and \A.3\ , 
and suppose that \9\ + 1/to ^ K and 1/ti ^ K\. Then with K = min(K, K\) we 
have, for any r ^ 2 and 5 > 0, 

\\^ r( l\\L 2 (dn) + II V<?IIl2(q) 

^c\\nv r q \\ L2{dn) + c(K,Vo\n) £ \\v s A q \\ L2{n) , (A.ie) 

s<r-l 

ll vr ^IL w + IIV r g|| L2(n) 

^5\\UV r q\\ L2{du) + C(l/6,K,Vom) IIVAgH^nj . (A.17) 

Lemma A. 7 (cf. [|, Proposition 5.9]). Assume that ^ r ^ 4. Suppose that 
\9\ ^ K and i\ ^ l/Ki, where i\ is as in Definition 3.5 of /fj/. If q = on dQ, then 
for m — 0,1, 

linv r g|| L2(5n) ^c(k, K 1 )(\\e\\ Loom + Yl P k ° 



k^r—2—m 

lL 2 (sn) • 



L 2 (dn) 

E H^IUnr ( A - 18 ) 



X 

fcsjr— 2+m 

//, in addition, \Vnq\ ^ s > and \Vnq\ ^ 2e || Vtv^zII L°°rsnv then 



V~ 2 e 



L 2 (dn) 



^ (A.19) 



:S2 



CHENGCHUN HAO AND T. LUO 



C K,K 1: - 



L°°(dn) 



fcsgr— 3 / k^r—1 



Lemma A. 8 (cf. ||, Proposition 5.10]). Assume that ^ r ^ 4 and that \6\ 
1/l ^ K. If q = on d£l, then 



L2(an) HWg|U (9n) + ||V^Ag|| L2(n) 



C (AT, Vol n, 



r-3 



x l|V^|| LOO(m) + ]T||V*Ag|| L , 



L2(fl) • 



(A.20) 



s=0 



Lemma A. 9 ( @, Lemma A.l]). If a is a (0, r) tensor, then with a = k/m and a 
constant C that only depends on m and n, such that 



V a 



if 



L s (dQ) 

m k m — k 



n II II I—a- II V7 m || a 

^ u 1 1 a 11^ (on) || v a \\ L p(dn) 



(A.21) 



2^p^s^q^ oo. 



s p q 

Lemma A. 10 ( 0, Lemma A.2]). Suppose that for i\ ^ 1/Ki 

\M{x\) —Af(x 2 )\ ^ e\, whenever \x\ — x~2\ ^ i\, Xx,X2 G 0*^, 

and 

where labiv) = 7ab(0, y). TTien if a is a (0, r) tensor, 

k 

* \\X7^\\ 

\Lp(dn) ' - ^ x- - ^, 



l a lll,(n-l)}>/(n-l-*p}(an) ^ ^(-^l) E ll^ a ll rpfflOt > 1 



n — 1 



£=0 



k-1 



l«llL-(sn) < 5 ll vfca Lp(an) + ^(#1) E ll V '"Ni. ■■„■-<>, 



a 



A; > 



n — 1 



V 



(A.22) 
(A.23) 

(A.24) 
(A.25) 



/or any 5 > 0. 



Lemma A. 11 ( Lemma A. 3]). With notation as in Lemmas and \A.1Q , we 

have 

k 



EH VJa L( n) < C 'll«lfc(n) EII V H 



LP(n) 



(A.26) 



Lemma A. 12 ( |jj Lemma A. 4]). Suppose that i\ ^ 1/K\ and a is a (0, r) tensor. 
Then 



|«|| L np/(n-fcp)(Q) ^C^^i 5 ^ ^ ||V 



£=0 

k 



lLp(n) 



£=0 



n 
p' 



(A.27) 
(A.28) 
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Lemma A. 13 ( @, Lemma A. 5]). Suppose that q = on dfl. Then 

ll9|lx 2(fi )^C(Vom) 1/n ||Vg|| i2(n) , (A.29) 

UVgll^n) ^C(Vom^ 2n \\Aq\\ L2{n) . (A.30) 
Lemma A. 14 ( |4], Lemma A. 7]). Let a be a (0,r) tensor. Assume that 

VolQ < V and \\9\\ Loc{m) + 1/lq < K, 
then there is a C = C(K, V,r,n) such that 

IMI L (n-iW(n-p) (an) ^ C||Va|| iP(n) +C\\a\\ LP{n) , l^p<n } (A.31) 

\\ v2a \\mn) < C (||nV 2 «|| i2(r _ 1}/ , I(an) + ||Aa|| £a(n) + ||Va|| x2(n) ) . (A.32) 
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